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Abstract 

We outline how unstable quantum fluctuations decohere classical fields in heavy ion collisions, 
leading to an equation of state and hydrodynamics. Explicit numerical realization of this frame- 
work in a scalar theory demonstrates that anomalously low values of ri/s can be generated. 



The motivation for this work is to achieve a deeper understanding of what one might con- 
sider the unreasonable effectiveness of hydrodynamics in heavy ion collisions. Hydrodynamics 
is a good effective field theory for the late time long wavelength behavior of a quantum field 
theory Why does it appear to work so well at times < 1 fm in heavy ion collisions ? An- 
other important motivation is to compute the right non-equilibrium initial conditions that can 
be matched on to viscous hydrodynamics, thereby eliminating an ad hoc feature of current phe- 
nomenology. 

Our approach is a weak coupling one, albeit the dynamics is very non-perturbative. Abun- 
dant analogies to such dynamics exist in other fields of physics. In this Color Glass Condensate 
eff'ective field theory f2], quantum fluctuations can be isolated and computed in principle order 
by order in as ; each order in this expansion includes resummations of different contributions 
depending on the nature of these fluctuations. Before the collision, we have to worry about fac- 
torization of quantum fluctuations into the wavefunctions of the incoming nuclei. The proper 
treatment of these is important to understand the energy evolution of the wavefunctions. Factor- 
ization implies that the rapidity F-dependence of the density matrix describing the n-body 
correlations of color charge densities p in the nuclear wavefunctions-described by the JIMWLK 
Hamiltonian |3|-is universal, regardless of whether the high energy probe of the nucleus is an 
electron, a proton or another nucleus. If true, this universality (proven currently ||4l only for 
"leading logs" in as Y at each order of perturbation theory) would be a powerful predictor of a 
wide range of phenomena in high energy QCD. 

The focus here is on heavy ion colUsions, where the leading order description is in terms of 
collisions of classical fields which best describe the high occupancy gluon fields in the nuclear 
wavefunctions |5]. In QCD one has classical gluon production, which gives the dominant con- 
tribution at early times Q. At early times, this matter, called the Glasma Q, contains lumpy 
(of size l/Qs, where Qs is the saturation scale) configurations of strong longitudinal chomo- 
electric and chromo-magnetic fields, giving rise to very anisotropic configurations with pressures 
Pt » /"l ~ IjS]. In gauge theories, these are known to give rise to instabilities, either of the 
Weibel f9] or Nielsen-Olesen type (10]. Therefore small 0(1) quantum fluctuations can grow 
to be as large as the background classical fields 1,11 J on parametric time scales (for expanding 
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systems) t ~ liv^il /as)/Qs ■ Since such contributions can occur at each order of perturbation 
theory, a resummation of these is required to achieve stable results. 

We showed recently that leading temporal instabilities (those that grow as (as exp(2 VGs t))", 
where n is an integer denoting the order in perturbation theory) can be resummed and expressed 
in terms of a gauge invariant spectrum of fluctuations on an initial Cauchy surface at t = 0^ [12l- 
This work suggests that the temporal evolution of inclusive quantities (such as the stress-energy 
tensor T^^ or correlators thereofi) can be expressed in terms of a "Master formula" 

<r''''>LLx+Li„st. = J [DpiDp2] W,, [pi] W,,[p2] J P'^] T^^mipuPi] + a](x) , (1) 

where the argument = (A, E) denotes collectively the components of the classical fields and 
their canonically conjugate momenta on the initial proper time surface; analytical expressions 
for these are available at t = 0^ |6|. Their temporal evolution is obtained by solving Yang- 
Mills equations [8 1. The Ws are the functional density matrices defined previously that obey the 
JIMWLK equation. 

The initial spectrum of fluctuations /^o^aj, Gaussian in the quantum fluctuations a, has a 
variance given by the small fluctuation propagator in the Glasma background field at t — > 0^. 
In practice, the path integral in a is determined by solving the classical Yang-Mills equations 
repeatedly with the initial conditions at t = 0^ given by 



init. init. 



d'ix) + c* d'^ix) 



(2) 



Here A'' denotes the quantum fields and their canonical conjugate momenta. The coefficients 
c^, with K denoting the quantum numbers labeling the basis of solutions, are random Gaussian- 
distributed complex numbers. Explicit expressions for the small fluctuations and their conjugate 
momenta, denoted here by a^(jc) were obtained in [12]. These need to be evaluated numerically- 
a challenging enterprise-but significant progress has been made in this direction. 

Considerable insight is gained by formulating the analogous problem in a massless scalar 
theory |13|. This theory, like classical Yang-Mills, is conformal and has instabilities due 
to parametric resonance of quantum fluctuations with the classical background. In the simple 
example of a fixed box, one sees that the leading order classical energy density and pressure 
for this theory are not single valued, with the pressure fluctuating rapidly as function of time. 
However, after adding fluctuations with Gaussian random coefficients (as in Eq. Q), the energy 
density and pressure rapidly develop a single valued relationship, namely, an equation of state 
(EOS). 

The EOS develops as a result of decoherence. Each of the trajectories corresponding to an 
initialization of the scalar analog of Eq. (|2| has a slightly different amplitude. In a conformal 
theory, different amplitudes have slighly different phases; for spatially independent fields and 
fluctuations in the theory, one can show that Tperiod - 18.2/§A0n,ax- In this case, the phase 
space density rapidly fills the constant energy Poincare surface uniformly; a simple exercise 
shows that this makes the stress-energy tensor traceless. (The canonical definition of the stress- 
energy tensor for the scalar theory is not traceless.) For an expanding l-nl-D scalar theory, one 
obtains Bjorken hydrodynamics. 

It is also instructive to look at spectral functions lfT4ll obtained from a Fourier transform of the 
imaginary part of the retarded Green functiorQ At early times, no quasi-particle behavior is seen. 



'The resummed expression for the latter is obtained from an equation analogous to Eq. jlj. The leading order 
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Figure 1 : Left figure: Time evolution of absolute values of diagonal elements the stress-energy tensor and the trace of the 
pressure tensor. Right figure: Evolution of the numerically extracted rjls ratio, compared to the perturbative value of this 
ratio, the value extracted from matching to first order viscous hydrodynamics and the conjectured universal lower bound 
rils = I /An. 



but it develops and one can extract a plasmon mass from the spectral function. Likewise, one 
obtains the occupation number from the resummed symmetric Green function (the sum of the 
Wightman functions G + + G+_). This develops a thermal structure = T jiuk - fj) - I /2, with 
the -1/2 denoting the contribution from vacuum fluctuations. Notably, there is a zero mode that 
is overoccupied relative to the thermal spectrum: it has been checked | fT4l that it demonstrates 
the characteristics of a Bose-Einstein Condensate, as also argued elsewhere 1 15 1. It is speculated 
that a transient condensate forms in the gauge theory analog as well 1 16 1. 

The (p* analogy to heavy ion collisions-of longitudinally expanding scalar fields-was ex- 
plored in IITtI . The theoretical framework is identical. One observes the pattern: Decoherence 
— » EOS — > Isotropization. Two striking features are illustrated in Fig. [T] Firstly, one notices 
(left figure) that despite a rapid red shift of the longitudinal pressure, the explosive instability 
growth allows the system to beat the expansion, leading to near isotropization at late times. It is 
important to note that the x-axis is in arbitrary lattice unit^ 

The other striking result is shown in Fig. [T] (right), where an effective extracted time de- 
pendent j]/s is shown. This is obtained by matching values of the diagonal components of the 
stress-energy tensor to first order viscous hydrodynamics using 

e 277 e 477 

P = - + — ■ P (3) 

^ 3 3t ' ^ 3 3t • ^ ' 

To extract rj/s, we use the Stefan-Boltzmann formula to estimate the energy and entropy density, 
e = and s = ^^^^ with s ~ e^^'*. In the hydrodynamical regime, 

(4) 

hydro 




expression is the fluctuation field that obeys the equation [n^ + V"(<f>{x'))] a(x) = 0. 

^Because of a logarithmic ultraviolet divergence in dynamically generated m^c/i^ terms in components of T''", com- 
posite operator renormalization is needed to relate computations at different lattice spacing and give physical meaning to 
time scales. Gauge invariance ensures this problem does not exist for QCD. 
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87 where A is the coefficient in the asymptotic behavior of the energy density, e ^ At ' . From 

88 this formula, we extract rj/s x 0.26. There is a large systematic uncertainty because (P^ - PJje 

89 at early times is best fit by a form exp(-Const.T^)-a faster relaxation of the pressure anisotropy 

90 than achieved through the power law behavior oc l/r^^^ typical in hydrodynamics. Nevertheless, 

91 the value extracted is two orders of magnitude lower the perturbative value [18] rjjs ~ lO^/g"* ~ 

92 40 for ^ = 4 and about three times the conjectured AdS/CFT bound. Such anomalously low 

93 viscosities also occur in turbulent phenomena [T9l where momentum transport occurs as if the 

94 viscosity were much smaller than transport cross-section estimates. 

95 It may be interesting to relate these results to work in the AdS/CFT framework describing 

96 potentially universal features in the relaxation of off-equilibrium strongly coupled systems to 

97 hydrodynamics li20J . Further, numerical simulations underway in the gauge theory case will 

98 improve state-of-the art Glasma-i-hydro computations 11211 by matching of the dynamics of the 

99 Glasma to viscous hydrodynamics. 
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